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In the first part of this review paper, the time-dependent Ginzburg-Landau theory is derived starting from 
the microscopic BCS model with the help of a derivative expansion. Special attention is paid to two space 
dimensions, where the entire crossover from the weak-coupling BCS limit to the strong-coupling BEC limit 
of tightly bound fermion pairs is accessible analytically. The second part deals with the dual approach to the 
. . . time-independent Ginzburg-Landau theory in three space dimensions. In this approach, the magnetic vortices 

0^ ' of a superconductor play the central role, and the superconductor-to-normal phase transition is understood as a 

proliferation of these vortices. 

I ! 

J:^ ! I. INTRODUCTION 

< ■ 

^ ' The role of topological defects in symmetry-breaking phase transitions cannot be underestimated. The prime example of a 
t defect-driven (equilibrium) transition is the superfluid-to-normal phase transition in a ''^He film, which can be understood as the 

unbinding of vortex-antivortex pairs [|l],^. Since the very existence of an ordered state and the presence of topological defects 
are both manifestations of spontaneously broken symmetries, it doesn't come as a surprise that various phase transitions can 
be understand in terms of defects. To appreciate this point of view, note that defects are regions, usually macroscopic in size, 
where the symmetry is broken differently from that in the rest of the system. Often, the symmetry is completely restored in the 
defect cores, meaning that they are in the normal state. (In systems like superfluid ■^He with a large symmetry group, the defect 
cores may be ordered states themselves, but of a different character.) Now, when defects with a normal-state core proliferate, the 
ordered state gets disordered and converted to the normal state. 

In recent years also the dynamics of phase transitions and the accompanying formation of defects have received much atten- 
tion (see Ref. for reviews). An important role is played here by the so-called time-dependent Ginzburg-Landau theory, which 
provides a phenomenological approach. Ideally, one would like to start with a microscopic model of the system under consid- 
eration and derive the time-dependent Ginzburg-Landau theory by integrating out irrelevant degrees of freedom, thus turning 
5^ the phenomenological theory into an effective one. For a superconductor this program, at least in the time-independent case, 
has been carried out by Gorkov [^. Starting from the microscopic theory of Bardeen, Cooper, and Schrieffer (BCS) [^], he 
^ derived by means of a Green function method the Ginzburg-Landau theory of superconductivity which had been proposed as a 
O phenomenological theory seven years before the BCS theory was formulated [||]. 

In the first part of this review paper, we reformulate Gorkov's derivation (extended to include time dependence) in terms of 
the functional-integral approach to quantum field theory. We will not only study the weak-coupling BCS limit of loosely bound 
Cooper pairs, but consider arbitrary values of the coupling constant Ao, including the limit Ao — oo, where the fermions 
^ form tightly bound pairs. When going from the weak-coupling to the large-Ag limit, the dynamics changes from dissipative to 
non-dissipative behavior 

In the second part of this review paper, we discuss the dual approach to the time-independent Ginzburg-Landau theory. It is 
a formulation directly in terms of magnetic vortices, in which the superconductor-to-normal is described as a proliferation of 
vortices. 



A. Notation 



We adopt Feynman's notation and denote a spacetime point by a; = = (t, x), jj. = 0,1, ■■■ ,d, with d the number of space 
dimensions, while the energy ko and momentum k will be denoted hy k — k^^ — (fco, k). The time derivative Oq = d/dt and the 
gradient V are sometimes combined in a single vector 9^ = {do, —V). The tilde on 9^ is to alert the reader for the minus sign 
appearing in the spatial components of this vector. We define the scalar product k ■ x — k^x^ = k^g^^k^, — k^t — k • x, with 
giiv — diag(l, • • • 7 — 1) and use Einstein's summation convention. Because of the minus sign in the definition of the vector 
9^, it follows that d^^a^ = 9oao + V • a, with an arbitrary vector 

Integrals over spacetime are denoted by 




Atd^x, 



1 



while those over energy and momentum by 



"27r (27rp" 



When no integration hmits are indicated, the integrals are assumed to be over all possible values of the integration variables. 
Natural units ?i = /cb = 1 are adopted throughout. 



II. TIME-DEPENDENT GINZBURG-LANDAU THEORY 



In this section, we derive the time-dependent Ginzburg-Landau theory for a superconductor, starting from the microscopic 
BCS model. Following Ref. [0], we integrate out the fermionic degrees of freedom in favor of a bosonic field describing 
the fermion pairs and use a derivative expansion method. We wish to obtain the Ginzburg-Landau theory not only in the weak- 
coupling BCS limit of loosely bound and overlapping Cooper pairs, but for arbitrary values of the coupling constant, in particular 
in the limit Ao — oo where the fermions form tightly bound pairs [^^. This crossover between the two limits has recently 
been studied in detail [[l0|-pj||. Since in two space dimensions this region becomes accessible by analytical methods, we shall, 
unless stated otherwise, restrict ourselves to d = 2 when carrying out momentum integrals. 



A. BCS Model 



The BCS Lagrangian reads [|l§] 



= : Cq + £i, 



(1) 



where £j = — AqV'i ip* V'i V'T is ^ contact interaction term, representing the effective attraction between electrons with bare 
coupling constant Aq < 0, and Co is the remainder. In Eq. (|l|), the field V'T(i) anticommuting field describing the electrons 
with mass m and spin up (down); zV) = e(— iV) — fiQ, with e(— iV) = — V^/2m, is the kinetic energy operator with the 
bare chemical potential /io subtracted. 

For computational convenience we introduce Nambu's notation and rewrite the Lagrangian (jl]) in terms of a two-component 
field 



(2) 



In order to integrate out the fermionic degrees of freedom, the zero-temperature partition function represented as a functional 
integral 



DV'TDV'exp [i £ , 



(3) 



must be written in a form bihnear in the electron fields. This is achieved by rewriting the quartic interaction term as a functional 
integral over auxiliary fields A and A*: 



(4) 



exp \^-iXo J i^i "01 ^1 
DA*DAexp 



The field equation for A*, 



I I ( A*ViV'T+V'T*V'IA-i^|A|^ 



A = XoiPiiPh 



(5) 



shows that the auxiliary field describes electron pairs. We will therefore refer to it as pair field. The partition function thus 
becomes 
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Z = I D7/>tD^ I DA*DA exp 



(6) 



X exp 



-A* ido + e(-*V) ) ^_ • 
Changing the order of integration and performing the Gaussian integral over the Grassmann fields, we obtain 

Z = JdA*BA exp (^*5eff[A*,A] + ^^|A|2^ , 

with Scs the one-loop effective action which, using the identity Det(yl) = exp [Tr In (A)], can be rewritten as 

5,ff[A*,A] = -zTrln(' P'> ' ^i^^ 7^ , 

where po — "ido and ^(p) e(p) — fiQ, with e(p) ~ p^/2m. 

The trace Tr appearing here needs some explanation. Explicitly, it is defined as 



(7) 



(8) 



S'cff — — iTr In [K{p, x) 



-itr In 



K{p,x)S{x - y)\ 



(9) 



where the trace tr is the usual one over discrete indices. We abbreviated the matrix appearing in (||) by K{p, x) so as to cover 
the entire class of actions of the form 
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ip^ {x)K(p, x)'ijj{x). 



The delta function in arises because K{p^ x) is obtained as a second functional derivative of the action 

5^S 



(5'0t(x) S^j;{x) 



K{p,x) S{x~y)\ 



(10) 



(11) 



each of which gives a delta function. Since the action has only one integral over spacetime, one delta function remains. 
Because it is diagonal, it may be taken out of the logarithm and the effective action (m can be written as 



S'cff = -itr / In [K{p, x)] 5{x - y) 



-itr 



y=x 



In the last step, we used the integral representation of the delta function: 



6{x) 



(12) 



(13) 



shifted the exponential function exp(zfc • y) to the left, which is justified because the derivative p^ does not operate on it, and, 
finally, set i/^ equal to . We thus see that the trace Tr in (^) stands for the trace over discrete indices as well as the integration 
over spacetime and over energy and momentum. The integral arises because the effective action (^ is a one-loop result with 
the loop energy and momentum. 

In the mean-field approximation, the functional integral (Q) is approximated by the saddle point: 



Z = exp iS'ofr[A*jf , A^f] 



-^0 Jx 



lA^fp 



where Amf is the solution of the mean-field equation 
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(14) 



(15) 



5A*{x) Ao 

If we assume the system to be spacetime independent, so that A^f (a;) — A, Eq. (Ol) yields the celebrated BCS gap equation 
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1 ^ /■ 1 



1 f 1 



2 A mr 

where ?/ is an infinitesimal positive constant that is to be set to zero at the end of the calculation, and 



(16) 



E{k) = ^e{^) + \^\^ (17) 

is the spectrum of the elementary fermionic excitations. 

B. BCS to BEC 

For a constant pair field, the effective action can be calculated in closed form. Writing 

f ko- e(k) -A A / fco - C(k) 

-A* fco + m y V ^0 + ^c^) 

and expanding the second logarithm in a Taylor series, we recognize the form 





( ' 




)- 


[ A* 





(18) 



5eff[A*,A] = 

where we ignored an irrelevant constant. The integral over the loop energy ko can be carried out to yield for the effective 
Lagrangian 



= ^ [E{k) - ^(k)] . 



(20) 



To this one-loop result we have to add the tree term |Ap/Ao. Expanding -E(k) in a Taylor series, we see that the effective 
Lagrangian also contains a term quadratic in A. This term diverges in the ultraviolet. To render the theory finite to this order, 
we have to introduce a renormalized coupling constant A defined by: 

' ' ' ^ ' (2.) 



To this order in the loop expansion there is no renormalization of the chemical potential, so that we can write jj, = fiQ. We 
regularize the diverging integral in Eq. ( pH ) by introducing a momentum cutoff A. In, for example, c? = 3, we then obtain 

where we omitted the (irrelevant) finite part of the integral. It should be remembered that the bare coupling constant Ao is 
negative, so that the interaction between the fermions is attractive. We can distinguish two Umits. One, the famous weak- 
coupling BCS limit, which is obtained by taking the bare coupling constant to zero, Aq 0~. Second, the limit which is 
obtained by letting Aq — oo. In this limit, the two-particle interaction is such that the fermions form tightly bound pairs of 
mass 2m 

To explicate this so-called Bose-Einstein condensation (BEC) limit in d = 2, we swap the bare coupling constant for a more 
convenient parameter, namely the binding energy Eq of a fermion pair in vacuum [ p3[ | . Both parameters characterize the strength 
of the contact interaction. To see the connection between the two, let us consider the Schrodinger equation for the problem at 
hand. In reduced coordinates, it reads 



V 



2 



TO 



+ Ao(5(x) 



^(x) = -eaV(x), (23) 



where the reduced mass is m/2 and the delta-function potential, with Ao < 0, represents the attractive contact interaction d in 
the BCS Lagrangian (|l]). We stress that this is a two-particle problem in vacuum; it is not the famous Cooper problem of two 
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interacting fermions on top of a filled Fermi sea. The equation is most easily solved by Fourier transforming it. This yields the 
bound-state equation 



1 ./^^.i (0)l..(^i), 



(25) 



where 1/(0) = m/27r is the two-dimensional density of states (per spin degree of freedom), and eA = A? /2m. This equation 
allows us to replace the bare coupling constant Aq with the binding energy ta- When substituted in the gap equation Eq. dl^), 
the latter becomes 



1 1 



(26) 



By inspection, it is easily seen that this equation has a solution 

A ^ 0, Mo ^ -5^a, (27) 

with a negative chemical potential. This is the strong-coupling BEC limit. To appreciate the physical significance of the specific 
value found for the chemical potential in this limit, we note that the spectrum i?b (q) of the two-fermion bound state measured 
relative to the pair chemical potential 2/Lto reads 

E^{q) = -Ea + ^ - 2^lo■ (28) 

The negative value for /io found in (|7|) is precisely the condition for a Bose-Einstein condensation of the composite bosons in 
the q = state — whence the name BEC limit. 

Since there are two unknowns contained in the theory, viz. A and /i, a second equation is needed to determine these variables 
in the mean-field approximation [^. It is provided by the requirement that the average fermion number N, which is obtained by 
differentiating the effective action (||) with respect to /i 

N^^, (29) 
be fixed. If the system is spacetime independent, this reduces to 

n^-itr I G{k)T3, (30) 

Jk 

where n = N/V, with V the volume of the system, is the constant fermion number density, T3 is the diagonal Pauli matrix in 
Nambu space. 



^3 = ( J _J ) , (31) 



and G{k) is the Green function. 



G(fc)=( ^° ] (32) 



-A* fco+^(k) 

" kl - E^k) +iT^\ A* ko e-^'^"" - e(k) 

Here, 77 is an infinitesimal positive constant that is to be set to zero at the end of the calculation. The exponential functions in 
the diagonal elements of the Green function are an additional convergence factor needed in nonrelativistic theories [n™. If the 



integral over the loop energy ko in the particle number equation (30) is carried out, it takes the familiar form 



The two equations ( [T^ and ( |30| ) determine A and 
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C. Derivative Expansion 



When the pair field A is spacetime-dependent, the integrals in ^ cannot be evaluated in closed form because the logarithm 
contains energy-momentum operators and spacetime-dependent functions in a mixed order. To disentangle the integrals we have 
to resort to a derivative expansion [ p^ in which the logarithm is expanded in a Taylor series. Each term contains powers of the 
energy-momentum operator which acts on every spacetime-dependent function to its right. All these operators are shifted to 
the left by repeatedly applying the identity 



f{x)p^g{x) = {pf, - idf,)f{x)g{x), 



(34) 



where f{x) and g{x) are arbitrary functions of spacetime and the derivative 9^ = {do, —V) acts only on the next object to the 
right. One then integrates by parts, so that all the p^'s act to the left where only a factor exp(ifc • x) stands. Ignoring total 
derivatives and taking into account the minus signs that arise when integrating by parts, one sees that all occurrences of (an 
operator) are replaced with fc^ (an integration variable). The exponential function exp(— ifc • x) can at this stage be moved to the 
left where it is annihilated by the function cxp{ik ■ x). The energy-momentum integration can now in principle be carried out 
and the effective action be cast in the form of an integral over a local density Cqb'- 



ScS 



-cff ■ 



(35) 



This is in a nutshell how the derivative expansion works [|17[|. 

To apply it to the BCS model and derive the Ginzburg-Landau theory we use the following decomposition in i 



-A* Po + ap)^ ° 



l-Gn 



A 
A* 



(36) 



where Go is the Green function ( [32| ) with A = 0. Apart from an irrelevant constant, this leads to the expression for the effective 
action 



Soff = "i Tr In 



1-Go 



STrV i 



Go 



A 

A* 

A 

A* 



(37) 



e=i 



For £ — I, the trace over the 2x2 matrix immediately yields zero. In a similar fashion all terms with £ odd, give zero. For 
the quadratic term, we obtain 



S. 



(2) 



i Tr 



■ A* 



Po + C(p) Po - S.{p) 



(38) 



where we recall the definition of the derivative p^ as operating on everything that appears to the right. Applying the derivative 
expansion rules outlined above, we can cast the quadratic term in the effective action in the form 



a(2) _ 
'eff — 



iTr 



1 



1 



ko + C(k) ko-ido~ C(k + iV; 



A*A. 



(39) 



Usually, the quartic terms are included in the Ginzburg-Landau theory without derivatives. We then can treat pq ± ^(p) as a 
c-number and 



2 [fc§-C^(k)]2 



(40) 



We may truncate the series in ( P7| ) here, provided the pair field A is small as is the case in the vicinity of the phase transition. 

To include the temperature T in the theory, we adopt the imaginary-time approach to thermal field theory [|8}|l^]. Very 
briefly, it can be derived from the corresponding quantum field theory at zero temperature simply by going over to imaginary 
times, t —IT, and substituting 



(41) 
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where g is an arbitrary function, while LUn denote the Matsubara frequencies. 



tt(3 ^2n, (bosonic) 
7r/3~^(2n + 1), (fermionic) 

with n an integer and f3 — 1 /T. With these rules, the Minkowski action S goes over into 



Jx Jq Jx 



(42) 



(43) 



where the superscript E on the action at the right-hand side is to indicate that it pertains to Euclidean rather than to Minkowski 
spacetime. 

For the case at hand we obtain in this way the finite-temperature action: 



^off - Tr 



1 



where the trace at finite temperature reads explicitly 



1 



2K + e^(k)]^ 



Tr : 







/3 r jdi, ^ 



(27r)° 



(44) 



(45) 



Let us for the moment consider only the time-independent part of the effective action ( |44| ) and expand in gradients. The sums 
over the Matsubara frequencies are carried out with the help of the formulas 



1 X 

2j 



n 



1 



1 



1 d 






(f) 


1 d 1 




49^ I 


S) 


X 


f3 d 







(46) 



where X and Y abbreviate the functions 



X tanh(/3C/2), Y ^ 1/ cosh^{(3^/2). 



(47) 



As was first shown by Drechsler and Zwerger [|T^, most of the momentum integrals can be performed analytically for arbitrary 
values of the coupling constant in two space dimensions. In this way, they arrived at the time-independent Ginzburg-Landau 
theory describing the crossover from the weak-coupling BCS limit to the strong-coupling BEC limit: 



with the coefficients 



C,s^cA*V'A + a\A\^-^b\A\\ 



1,(0 \ /Ae''" 

a=^-\ni — ] e{fji) + ln(/3e,/4) + ln(/3|Ai|/2) tanh(^///2) 



+sgn(^) / dx \n{x)Y{x) 

Jp\^\/2 



(48) 



(49) 



b = 



YS) 1 r2 roc 

/320(^) - — tanh(/3^/2) + sgn(/.)^ / dx 



7C(3) 
2 



Xix) 



(50) 



and 



7 



8m 



7C(3) 
27r2 



fl2 rco 

4 i/3|M|/2 



X{x) 



(51) 



Here, we introduced the integration variable x — /3^/2, 7 = 0.577216 • • • is Euler's constant, 0{x) is the Heaviside unit step 
function, sgn(a;) the sign function, and C,{x) Riemann's zeta function, with C (3) = 1.20206 • • •. (The second term at the right- 
hand side of Eq. (pO| ) differs slightly from the corresponding term in Ref. ([l^], but is consistent with a later work by Zwerger 
and collaborator [[20|].) 

When studying the time-dependence of the effective action (p4|), special care has to be taken with analytic continuation. We 
follow Sa de Melo, Randeria, and Engelbrecht [ 12 1 and first analytic continue, using the formula 



1 



1 



1 



2^ + 90 



i7r(5(2^ + go) 



(52) 



before expanding in time derivatives. (In Ref. [ |lO[|, the expansion was done first, leading to results which are not consistent with 
those known in the BCS limit [ pl] , p2[ ].) In Eq. ([52|), P stands for the principal part, while ui is a bosonic Matsubara frequency, 
which at the right-hand side is analytic continued to the real axes by replacing iw; with go + *'7- In this way, we find for the 
dynamic part (in Minkowski spacetime) 



Cdyn = [Q'{id^) - i7rQ"(i9o)]A*A, 



(53) 



with 



3' (go) 
?"(go) 



90 



-X 



,k 2^(2^ + 90) 
i^/(0)tanh(/3go/4)^^(Ai-9o/2). 



(54) 



The most important result to be noted here is that the low-energy dynamics is dissipative when the chemical potential /i is 
positive [ 12 1. This is because the pairs can break up and decay into a continuum of fermionic excitations. On the other hand, 
for negative values of /i, where the pairs are more tightly bound, the time-dependent Ginzburg-Landau theory describes a purely 
propagating pair mode. To investigate this point further, let us consider the two limits in detail. 



D. BCS limit 



In the weak-coupling BCS limit, where the chemical potential is well approximated by the Fermi energy /i — n,p 
k-p the Fermi momentum, and eo/ep 0, we recover the standard result (adjusted for the reduced space dimensionality) 



fc§/2m, with 



£eff = J^(O) <^ A* 



Inl I I - 



8T0 



with VF — k-p/m the Fermi velocity, ^0 the BCS correlation length 



A-^|A|' 



^0^ 



7C(3) vl 



487r^ 



rp2 ' 

-'o 



and Tq the BCS transition temperature [ p^ ] 



To = — (2^ea) 

TT 



1/2 



(55) 



(56) 



(57) 



expressed in terms of the binding energy ta- Comparing the two terms in Eq. ( |55| ) involving time derivatives, we recognize a 
series expansion in powers of Qq/T; the terms without derivatives constitute an expansion in powers of | Aj/T [^. It therefore 
follows that Eq. ( p5| ) represents the first terms in a high-temperature expansion. Since we neglected time derivatives in the 
quartic and higher order terms, we are implicitly assuming that qo >> | A| [p2|]. The time-dependent Ginzburg-Landau theory 
(^) obtained in the weak-coupling BCS limit is of the dissipative type often used in the context of defect formation in a 
symmetry-breaking phase transition (for reviews see Ref. [^). For applications in the context of superconductivity see, for 
example, Refs. 
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E. BEC Limit 



In the strong-coupling BEC limit, the general form of the effective action ( |4q ) and ( p3| ) reduces to [ p^ 

£eff = A* 6ao + Mb + ^) A-Ab|A|^ (58) 



where we introduced a rescaled pair field 



.(0) 



4m 



1/2 



The effective theory (|8|) is precisely of the form of a Gross-Pitaevski theory [^t]], describing a weakly interacting composite 
Bose gas with a mass 2m (as expected), a small chemical potential which vanishes on approaching the critical temperature from 
below [see Eq. (|7|)] 

Aib = 2\fi\ In (^) , (60) 



or using Eq. 



and a repulsive contact interaction 



2|mI 



lAP 

Mb = — , (61) 



independent of the binding energy characterizing the interaction between the electrons. This is special to two dimensions; in 
d space dimensions we have instead 

1 J /A l-d/2 

^ r(2 - d/2) m^/2 ' ^^^^ 



with r(a:) the Gamma function. Note the absence of any temperature-dependence in the effective theory (p8|). At zero temper- 
ature, exactly the same effective theory was obtained by Hausmann, using a self-consistent Green function method JTl]]. (See 
Ref. for a derivation along the lines presented here.) This is, as that author argued, because in the BEC limit, the critical 
temperature Tq is much smaller than the dissociation temperature Tdiss ~ Ca at which the tightly bound fermion pairs are broken 
up by thermal fluctuations. Hence, for all temperatures in the range T <Tq « Tdiss we are effectively in the zero-temperature 
regime Jill]. 

The Gross-Pitaevski theory ( ^8| ) describes a gapless, purely propagating mode, viz. the Goldstone mode associated with the 
spontaneously broken global U(l) symmetry — the so-called Anderson-Bogoliubov mode [^. 

The spacetime-dependent effective theory of a superconductor can only be derived in a few special cases [ p2[ ] : one being close 
to the transition temperature where the expansion parameter is 1 /T and another being close to the absolute zero of temperature 
where the expansion parameter is 1/|A|. Outside these regimes, the effective theory depends on the ratio |V|/9o, and an 
expansion in both time derivatives and gradients is not possible. It is amusing to note that in going from the weak-coupling BCS 
limit to the strong-coupling BEC hmit, we move from one valid regime to the other. 



III. DUAL THEORY 



In this section we investigate the dual formulation of the time-independent Ginzburg-Landau theory. This approach, in which 
the magnetic vortices of a superconductor play the central role, originates from lattice studies that started more than two decades 
ago [p9[-^6|. These were in turn instigated by the success of the Kosterlitz-Thouless theory describing the phase transition in 
a superfluid film as the unbinding of vortex-antivortex pairs [|l]J^]. In the dual formulation of the Ginzburg-Landau theory, the 
superconductor-to-normal phase transition is understood as a proliferation of magnetic vortices. A detailed presentation of these 
matters as well as an extensive list of references to the literature can be found in Ref. [^. 
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A. Electric Current Loops 



To account for the magnetic interaction we couple the Ginzburg-Landau theory in the usual, minimal way to a vector potential 
A. We also rescale the pair field such that /3 times the Hamiltonian becomes 



Wgl- |(V-2ieA) 



1 



(V X A) 



1 

2^ 



(V • Af 



(64) 



where we added a gauge-fixing term with parameter a. To acquire a physical understanding of what this Hamiltonian describes 
p^], we recall that a |^|^-theory gives a field theoretic description of strings with contact repulsion p^. This equivalence rests 
on Feynman's observation [Eop that the Green function 



G(x) 



k2 



(65) 



can be expressed as a path integral. An easy way to see this is to invoke Schwinger's proper-time method, which is based on 
Euler's form 



1 



1 



At 



to write the right-hand side of (^5|) as [41 1: 

G(x 



nz) Jo 



dr e 



dr e 



T e 



gik-Xg-rk 



(66) 



4ttt J 



3/2 



-ix^ /t 

e 4 ' . 



(67) 



The factor (I/Attt)^/^ exp(— ix^/r) appearing here can be interpreted as describing a Brownian string trajectory, showing that 
one endpoint of the string (located at x) has a Gaussian distribution with respect to its other endpoint (located at the origin) [|4^]. 
If we imagine the string to be composed of N links, each of length a, then in the limit where N oo, a ^ 0, the variable r 
parameterizing the string stands for r = No? The integration over r in ( ^7[ ) indicates that the string can be arbitrary long, 
but the weighing factor exp(— m^r) exponentially suppresses long ones in the normal phase where > 0. To understand the 
physical meaning of this factor, let us return to the discrete string model and write 



so that it becomes 



= 6[o-(a) - o-cr]/c 



(68) 



(69) 



with (Jeff (a) = a{a) — cjcr the effective string tension, and L — Na the length of the string. The continuum limit is obtained by 
simultaneously letting a ^ and (T(a) o'er, in such a way that the right-hand side of (68) tends to the finite value m|. We 
thus see that the factor cxp(— m^r) weighs strings according to their lengths. 

Following Feynman [|40|], we write the right-hand side of (^) as a path-integral, i.e., as a sum over all possible string trajec- 
tories having one endpoint at x(0) — and the other at x(r) ~ x [p3|]: 



G(x) 



/'x(r)— X 

dr / Dx(T')e-^", 

Jx(0)=0 



with the (Euclidean) action 



So^ rdr'[ix2(r')+m^], 
Jo 



(70) 



(71) 



where x(t) = dx(T)/dT. 

In a similar way, also the partition function of the free theory, which written as a functional integral reads 



Zo 



D(t)*D(j) exp 



Ho 



(72) 
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can be be represented as a path integral, this time involving only closed strings: 



In(Zo) = - ln[Dct(p2 + m^)] = -Trln(p2 + m^) (73) 

-So 



Here, we used Euler's form ( |66| ) in the limit of small z. An additional factor 1 /r arises because one can start traversing a closed 
string anywhere along the loop. The |0|^ -interaction in the Ginzburg-Landau model can be shown to result in the additional term 
0] 

Sx = -A f dr/dr^ 5 [x(r/) - x^)] (74) 
Jo 

in the action, which gives an extra weight each time two strings — one parameterized by r/ and one by — intersect. Physically, 
it represents a repulsive contact interaction between strings. Finally, the coupling of the field (j) to the magnetic vector potential 
A via the electric current 

je = -2ei(0*V(/) - (j)V<j)*) - 2(2e)2A|(/)p (75) 
with a charge 2e, results in the extra term |j3^] 

Se^2ie [ dT'i(r') • A[x(t')], (76) 

showing that the strings described by the Ginzburg-Landau theory carry an electric current. Pasting the pieces together, we 
conclude that the partition function of the Ginzburg-Landau theory can be equivalently represented as a grand canonical ensemble 
of fluctuating electric current loops, of arbitrary length and shape [p^]: 



N=0 ■ 1=1 



e~""^^ (77) 



with the action 



N n 

Sgl^Y. dr/ (ix^W) + + 2zex(r/) • A[x(r/)]} 
1=1-^° 

+X f2 r dr/ r d4 S [x(r/) - x(r^)] . (78) 



Lk=l 



On entering the superconducting phase, characterized by a sign change in the mass term of the Ginzburg-Landau theory, the 
effective string tension approaches zero as 



aeff(a)cxm^cx(r-re)"", (79) 

with h' the correlation length exponent, and the electric current loops proliferate. This proliferation is also signaled by the 
absolute value |(/)| of the pair field which then develops a vacuum expectation value, identifying it as the order parameter 



B. Disorder Parameter 

The dual formulation of the Ginzburg-Landau gives a similar representation of the partition function as the path integral (|77|); 
however this time not in terms of electric current, but of magnetic vortex loops. To derive it, we start with a trick [ p^ and 
introduce a (hypothetical) magnetic monopole at some point z inside the superconductor. A monopole is a source of magnetic 
flux. Due to the Meissner effect, the flux lines emanating from the monopole are squeezed into a flux tube. In this way, we 
managed to create a magnetic vortex at zero external field. Electrodynamics in the presence of a monopole was first described 
by Dirac | p7[ | who argued that the physical local magnetic induction h is given by the combination V x A(x) — B^(x). The 
subtracted plastic field 
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(80) 



with ^0 ~ n/e the magnetic flux quantum, removes the field of the so-called Dirac string running along some path from the 
location z of the monopole to infinity. On account of Stokes' theorem, the plastic field satisfies the equation 



V-BP(x) = -p,n(x), (81) 



with /Om(x) — ^0 d{x — z) the monopole density. 
We continue by writing (f>{x.) as 



^)^^g2^v(x)^ (82) 



where is a constant solution of the mean-field equation. This approximation, where the phase of the order field is allowed to 
vary in space while the modulus is kept fixed, is called the London limit. The Ginzburg-Landau Hamiltonian then becomes after 
integrating out the phase field (p: 

n^GL - ^(V X A - BP)2 + ImiA. - A, + ^(V . A)^ (83) 



where 



TUA = 2e|(/)| (84) 



is the inverse magnetic penetration depth ttia = • We gave H the superscript P to indicate the presence of the monopole. 
While the Dirac string is immaterial in the normal phase, Nambu [ p^ argued that it acquires physical relevance in the supercon- 
ducting phase where it serves as the core of the magnetic vortex or Abrikosov flux tube originating at the monopole. 

The energy Ey of this configuration is easily obtained by substituting the solution of the field equation for the gauge field 



A,(x)= / G,,(x-y) [VxBP(y)]^., 
•'y 



(85) 



with Gij the gauge-field Green function 



back into the Hamiltonian (|3]). The energy is divergent in the ultraviolet because in the London limit, where the mass jm^j of 
the (/)-field is taken to be infinite, the vortices are considered to be ideal lines. For a finite mass, a vortex core has a typical width 
of the order of the coherence length ^ — l/\m^\. This mass therefore provides a natural ultraviolet cutoff to the theory. Omitting 
the irrelevant (diverging) monopole self-interaction, one finds [ p^ ] 

Ev = ^.g' ^ da;, ^ dy., G(x - y) = ayL,. (87) 

Here, 

G(x) = / — 2- = n— (88) 

Ak^+m^ 47r |x| 
is the scalar Green function, g abbreviates the combination 

9 = $oTOA, (89) 
while Lz denotes the (infinite) length of the Dirac string and 

its line tension, first calculated by Abrikosov [S^. The value kql = 1/ of the Ginzburg-Landau parameter 
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■niA 



Al 



or 



2^' 



(91) 



separates the type-II regime (kql > l/\/2), where isolated vortices can exist, from the type-I regime (kql < l/\/2), where 
a partial penetration of an external field is impossible. Remembering that was the Dirac string, we see from ( p7| ) that in the 
superconducting phase it indeed becomes the core of the magnetic vortex originating at the monopole, as was first observed by 
Nambu Q. 

The operator T^(iz) describing the monopole with its emerging flux tube 51 1 is easily constructed by noting that in 
the functional-integral approach, a given field configuration is weighted with a Boltzmann factor exp (— J^^gl)' ^^^^ the 
Hamiltonian given by Eq. (|83|). From this we infer that the explicit form of the vortex operator is 



V{L^) = exp J / (V X A) • - i (B^ 



(92) 



We next demonstrate that this operator can be used to distinguish the superconducting from the normal phase [ pO|pl| ]. To this 
end, let us consider the correlation function {V {Ly)V* {L^)) , where F*(Lz) describes an additional antimonopole brought into 
the system at z, with being the accompanying Dirac string running from infinity to z. Since all the integrals involved are 
Gaussian, this expectation value can be evaluated directly. We proceed, however, in an indirect way for reasons that will become 
clear when we proceed, and first linearize the functional integral over the gauge field by introducing an auxiliary field h. In the 
gauge V • A = 0, which corresponds to setting a = in the Hamiltonian (p3|), we find [p2|] 



' 2 



{V{L^)V*{L-^)) = 

DADh exp -ih2 + ih (VxA 

where the plastic field satisfies ( pTj ) with the monopole density given by 

Pm(x) = $o['5(x - z) - <5(x - Z)]. 

To appreciate the physical relevance of the auxiliary field, let us consider its field equation 

h = i(V X A-RP) = ih.. 

It tells us that apart from a factor i, the auxiliary field h can be thought of as representing the local magnetic induction h. 
The integral over the vector potential is easily carried out by substituting the field equation for A, 

A = X h, 



(93) 



(94) 



(95) 



(96) 



back into (^3|), with the result 

{V{L^)V*{L^)) 



Dh exp 



-(V X h)2 + h2 



(97) 



This shows that the magnetic vortex described by the plastic field B^ couples to the fluctuating massive vector field h, with a 
coupling constant g introduced in Eq. (|9|). As the temperature approaches the critical temperature from below, (j) tends to zero, 
so that the vortex decouples from the auxiliary field h. The finite magnetic penetration depth in the superconducting phase is 
reflected by the mass term in (^tJ). 

After carrying out the integral over h in (p7|), we obtain for the correlation function 

= (98) 

exp|-i ^ [An(x) G(x - y) p„,(y) + m^Sf (x) G(x - y) i?f (y)] 

The first term in the argument of the exponential function contains a diverging monopole self-interaction for x = y. This 
divergence is irrelevant and can be eliminated by defining a renormalized operator 

K (iz) = V{L^) exp [i$gG(0)] . (99) 

The second term in the argument is the most important one for our purposes. It represents a Biot-Savart interaction between two 
line elements dxi and dyi of the magnetic vortex (see Fig. Rb. 
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FIG. 1. Biot-Savart interaction (wiggly line) between two line elements dxi and dyi of a magnetic vortex (straight line). 



For the renormalized operators we find 



{Vr{L^)V*{L^)) = cxp(-(Tyi2;z) exp 



<i>^ e 



47r 



(100) 



where L^z denotes the length of the flux tube connecting the monopole at z with the antimonopole at z. Initially, the two Dirac 
strings may run to infinity along two arbitrary paths. Due to the string tension, however, they join on the shortest path L^z 
connecting the monopoles. 



The coiTelation function (100) behaves differently in the two phases. In the superconducting phase, where rriA ^ 0, the 



first factor dominates. It represents the confining linear potential between the monopole and the antimonopole. As a result, the 
correlation function decays exponentially for distances larger than 1 /ay- 



(101) 



On the other hand, in the normal phase, where the gauge field is massless, the confinement factor in the correlation function 
(IOC) disappears, while the argument of the second exponential turns into a pure Coulomb potential. The coiTelation function 



remains, consequently, finite for large distances: 



{Vr{L,,)V;{L,,)) 



1, 



(102) 



thus indicating a proliferation of magnetic vortices in the normal phase. This demonstrates that the vortex operator can be used 
as an order parameter to distinguish the two phases. Since it develops a vacuum expectation value not in the superconducting, 
but in the normal phase it is refeiTed to as a disorder parameter 

It is interesting to consider in detail the limit T t Tc, where the coupling constant g = ^aniA tends to zero and the magnetic 
vortex decouples from the massive vector field [^]. In Eq. (97), this limit yields the constraint V x h = which can be solved 
by setting h = V7. The correlation function then takes the simple form 



{v,{z)v;{i)) = 



J D7exp -i J (V7)'+i J 



IPn 



(103) 



In the absence of monopoles, the theory reduces to that of a free gapless mode 7 that may be thought of as representing the 
magnetic scalar potential: 



V7 = iV X A. 

This follows from combinin g the physical interpretation of the vector field h ( p5| ) with the equation h = V7. 
The correlation function ( 103 ) can be put in the form 



(K(z)v;.*(z) 



,i$o[7(z)-7(z)] 



(104) 



(105) 



where the average at the right-hand side is taken with respect to the free scalar theory. This equation shows that in the normal 
phase (T > Tc), the Dirac string looses its physical relevance, the right-hand side depending only on the end points z and z, not 
on the path L^z connecting these points. The notion of a magnetic vortex is of no relevance in this phase because the vortices 
proliferate and carry no energy. This is the reason for omitting any reference to vortex lines in the argument of V in Eqs. (103) 
and dlO^). 



C. Magnetic Vortex Loops 

From the above results for a single vortex, we can now easily infer the dual formulation of the Ginzburg-Landau theory. In 
this formulation, the partition function is written as a grand canonical ensemble of fluctuating magnetic vortex loops, of arbitrary 
length and shape: 
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N=0 ' 1=1 



n 



Dx(r/) 



with the dual action, [cf. Eq. ([T^] 

^ rn 

^duai = E / "^^l {i^'W) + + **oi(r/) • h[x(T/)]} 
+uf] Tdr/ rdr^<5[x(r/)-xK)]. 



-Sdua 



(106) 



(107) 



We have included here a mass term representing the intrinsic vortex line tension [cf. Eq. (pq)], and also a contact repulsion 
between the vortices parameterized by u. 

The equivalent field representation of this dual theory reads [ ^3|j54| , p7| , ^5|j52| , F6| l. 



/ 



Z = / DhDV'*DV' exp - / Hd 



dual 



with the Hamiltonian 



ual 



-(V X h)2 



2ml '2 



|(V-i$oh)^Ap+m^|V'P + u|V'l' 



(108) 



(109) 



where the disorder field is minimally coupled to the massive vector field h, representing the local magnetic induction. The 
dual formulation contains the same information as the original, Ginzburg-Landau formulation. For example, the vortex line 
tension ( 90 ) appears in the dual theory as a one-loop on-sh ell rn ass correction stemming from the graph depicted in Fig. [l| which 
we now interpret as a Feynman graph of the dual theory (109), with the straight and wiggly line denoting respectively the i/;- 
and h-field Green function. Also the fixed points of the Ginzburg-Landau theory map onto those of the dual field theory [^]. 
And both formulations can be used to study the critical behavior of the superconductor-to-normal phase transition (see Refs. 
||3^,0|7|g|6|g,|8|l for the dual approach). 

Acknowledgments 

It is a pleasure to thank N. Antunes, L. Bettencourt, A. Leggett, D. Steer, and G. Volovik for useful discussions during the NATO 
Winter School and European Science Foundation (ESF) Workshop Topological Defects and the Non-Equilibrium Dynamics of 
Symmetry Breaking Phase Transitions, Les Houches, February 16-26, 1999, and M. Crisan for helpful correspondence. 

This work is performed as part of a scientific network supported by the ESF (see network's URL, http://www.physik.fu- 
berlin. de/^defect) . 



[1] V. L. Berezinskii, Sov. Phys. JETP 34, 610 (1972). 

[2] J. M. Kosterlitz and D. J. Thouless, J. Phys. C 6, 1181 (1973). 

[3] NATO Advanced Study Institute Formation and Interactions of Topological Defects edited by A.-C. Davis and R. Brandenberger, (Plenum 

Press, New York, 1995) 
[4] L. P Gorkov, Sov. Phys. JETP 9, 1364 (1959). 

[5] J. Bardeen, L. N. Cooper, J. R. Schrieffer, Phys. Rev. 108, 1175 (1957). 

[6] V. L. Ginzburg and L. D. Landau, Zh. Eksp. Teor. Fiz. 20 (1950) 1064; reprinted in L. D. Landau, Collected Papers (Pergamon, London, 
1965), p. 546. 

[7] A. M. J. Schakel, On broken symmetries in Fermi systems. Ph. D. Thesis, University of Amsterdam (1989). 
[8] D. M. Eagles, Phys. Rev. 186 (1969) 456. 

[9] A. J. Leggett in Modern Trends in the Theory of Condensed Matter, edited by A. Pekalski and J. Przystawa (Springer- Verlag, Berlin, 
1980), p. 13. 

[10] M. Drechsler and W. Zwerger, Ann. Phys. (Germany) 1, 15 (1992). 
[1 1] R. Haussmann, Z. Phys. B 91, 291 (1993). 

[12] C. A. R. Sa de Melo, M. Randeria, and J. R. Engelbrecht, Phys. Rev. Lett. 71, 3202 (1993). 



15 



M. Marini, F. Pistolesi, and G. C. Strinad, Eur. Phys. J. 1, 151 (1998). 

V. N. Popov, Functional Integrals and Collective Excitations (Cambridge University Press, Cambridge, 1987). 
M. Randeria, J. -M. Duan, and L. -Y. Shieh, Phys. Rev. B 41, 327 (1990). 

R. D. Mattuck, A Guide to Feynman Diagrams in the Many-Body Problem (McGraw-Hill, New York, 1976). 
C. M. Fraser, Z. Phys. C 28, 101 (1985); I. J. R. Aitchison and C. M. Fraser, Phys. Rev. D 31, 2605 (1985). 
R. J. Rivers, Path Integrals in Quantum Field Theory (Cambridge University Press, Cambridge, 1987). 
J. I. Kapusta, Finite-Temperature Fi eld Theory (Cambrid ge University Press, Cambridge, 1989). 



S. Stintzing and W. Zwerger, e-print :ond-mat/970312g (1997). 
A. Schmid, Phys. Kond. Materie 5, 302 (1966). 
E. Abrahams and T.Tsuneto, Phys. Rev. 152, 416 (1966). 
A. M. J. Schakel, Int. J. Mod. Phys. B 8, 2021 (1994). 

PG. de Gennes, Superconductivity in Metals and Alloys (Benjamin, New York, 1966). 

M. Tinkham, Introduction to Superconductivity (McGraw-Hill, New York, 1975). 

M. Crisan, Theory of Superconductivity (World Scientific, Singapore, 1989). 

E. P Gross, Nuovo Cimento 20, 454 (1961); L. P Pitaevskii, Sov Phys. JETP 13, 451 (1961). 

A. M. J. Schakel, in Correlations, Coherence, and Order, edited by D. V. Shopova and D. I. Uzunov (Plenum Press, New York, 1999), p. 
295. 

T. Banks, B. Meyerson, and J. Kogut, Nucl. Phys. B 129, 493 (1977). 

M. Peskin, Ann. Phys. 113, 122 (1978). 

P R. Thomas and M. Stone, Nucl. Phys. B 144, 513 (1978). 

W. Helfrich and W. Miiller, in Continuum Models of Discrete Systems (Waterloo University Press, Waterloo, 1980), p. 753. 
C. Dasgupta and B. I. Halperin, Phys. Rev Lett. 47, 1556 (1981). 
H. Kleinert, Lett. Nuovo Cimento 35, 405 (1982) 
J. Bartholomew, Phys. Rev. B 28 (1983) 5378. 
R. Savit, Rev. Mod. Phys. 52, 453 (1989). 

H. Kleinert, Gauge Fields in Condensed Matter (World Sc ientific, Singapore, 1989) , Vol. I. 



A. M. J. Schakel, Boulevard of Broken Symmetries, e-print cond-mat / 9805152 (1998). 

K. Symanzik, in Euclidean Quantum Field Theory, edited by R. lost (Academic, New York, 1969). 

R. P Feynman. Phys. Rev. 80, 440 (1950). 

J. Schwinger, Phys. Rev. 82, 664 (1951). 

H. Kleinert, Path Integrals in Quantum Mechanics, Statistics and Polymer Physics, 2nd Edition (World Scientific, Singapore, 1995). 
R. P Feynman, Rev. Mod. Phys. 20, 367 (1948). 

G. Parisi, Statistical Field Theory (Addison- Wesley, New-York, 1988). 

E. Copeland, D. Haws, S. Holbraad, and R. Rivers in The Formation and Evolution of Cosmic Strings, edited by G. W. Gibbons, S. W. 

Hawking, and T. Vachaspati (Cambridge University Press, Cambridge, 1990), p. 35; ibid p. 49. 

M. Kiometzis, H. Kleinert, and A. M. J. Schakel, Fortschr. Phys. 43, 697 (1995). 

P A. M. Dirac, Phys. Rev 74, 817 (1948). 

Y. Nambu, Phys. Rev. D 10, 4262 (1974). 

A. A. Abrikosov, Sov Phys. JETP 5, 1 174 (1957). 

E. C. Marino, Phys. Rev D 38, 3194 (1988); E. C. Marino, G. C. Marques, R. O. Ramos, and J. S. Ruiz, Phys. Rev. D 45, 3690 (1992); 
E. C. Marino, Ann. Phys. (NY) 224, 225 (1993). 

A. Kovner, B. Rosenstein, and D. Eliezer, Mod. Phys. Lett. A 5, 2733 (1990); Nucl Phys. B 350, 325 (1991); A. Kovner and B. Rosenstein, 
Phys. Rev. Lett. 67, 1490 (1991). 

M. Kiometzis and A. M. J. Schakel, Int. J. Mod. Phys. B 7, 4271 (1993). 
K. Bardakci and S. Samuel, Phys. Rev D 18, 2849 (1978). 

H. Kawai, Prog. Theor Phys. 65, 351 (1981). 

A. Kovner, P. Kurzepa, and B. Rosenstein, Mod. Phys. Lett. A 8, 1343 (1993). 
M. Kiometzis, H. Kleinert, and A. M. J. Schakel, Phys. Rev Lett. 73, 1975 (1994). 

I. F Herbut, J. Phys. A: Math. Gen. 29, 423 (1996). 



C. de Calan and F. S. Nogueira, e-print :ond-mat/ 990324 7 (1999) 



16 



